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We study the quantum dimer model on the triangular lattice, which is expected to describe the 
singlet dynamics of frustrated Heisenberg models in phases where valence bond conhgurations dom- 
inate their physics. We hnd, in contrast to the square lattice, that there is a truly short ranged 
resonating valence bond (RVB) phase with no gapless collective excitations and with deconhned, 
gapped, spinons for a finite range of parameters. We also establish the presence of three crystalline 
phases in this system. 



PACS numbers; 75.10.Jm, 74.20.Mn 

The search for an RVB phase of frustrated magnets, 
inspired by ideas of Pauling and begun in earnest by 
Anderson has been one of the recurrent themes in 
research on the cuprate superconductors over the past 
decade and somewhat more. Shortly after Anderson's 
1987 paper on the cuprates 1^, RVB theory bifurcated — 
with a gapless RVB state (with valence bonds on many 
length scales) pursued largely by means of gauge the- 
oretic treatments introduced by Baskaran and Ander- 
son Q defining one track, while Kivelson, Rokhsar and 
Sethna jsj hewed closely to the original vision and pur- 
sued the study of a short ranged RVB state. The latter 
proposal, of a state with exponentially decaying spin-spin 
correlations and no long range valence bond order, was 
expected to lead to a gapped collective excitation spec- 
trum and deconfined, gapped, spinons. Unfortunately, it 
turned out that this was hard to arrange on the square 
lattice — the simplest implementation of RVB ideas, the 
quantum dimer model exhibits crystalline order and 
confined spinons, except at a critical point This fact 
is of a piece with an instability of the paramagnetic phase 
of the 0(3) non-linear sigma model to breaking transla- 
tional symmetry due to Berry phase effects |^] . 

In this Letter, we report that the simplest quantum 
dimer model on the triangular lattice does possess a 
short ranged RVB phase with gapped collective modes, 
gapped deconfined spinons and spin-charge separation in 
its charged excitation spectrum. We establish the pres- 
ence of nearby crystalline phases with confined spinons. 
We also suggest that a connection, made previously by 
P. Chandra and ourselves between quantum dimer 
models and frustrated transverse field Ising models can, 
in this problem, be plausibly extended to track the tran- 
sition out of the RVB phase. This conjecture is closely 
connected to the spinon deconfinement mechanism pro- 
posed by Read and Sachdev 0,|ll|, Wen [|| and the 
ideas of Senthil and Fisher Q . 

We study the Rokhsar-Kivelson quantum dimer Hamil- 
tonian generalized to the triangular lattice (as the con- 
trasts are instructive, we will comment on the square 
lattice results along the way) : 



H = -tf + vV = J2\-tJ2[\ ^){\\\ + h.c.) (1) 
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Here, the sum on i runs over all of the Np plaquettes, and 
the sum on a over the three different orientations of the 
dimer plaquettes, namely V!i rotated by and ±60°. We 
refer to the plaquettes with a parallel pair of dimers as 
fiippable plaquettes. As a complete orthonormal basis set 
we use {|c) | c = l...A''c}, where |c) stands for one of the 
Nc possible hardcore dimer coverings of the triangular 
lattice. V is diagonal in this basis, with V\c) = nfi{c) \c) 
measuring the number, nfi(c), of fiippable plaquettes in 
configuration c. 

Rokhsar and Kivelson derived the square lattice ver- 
sion of H as the leading effective Hamiltonian in the sin- 
glet manifold consisting of nearest neighbor valence bond 
coverings of the lattice by utilizing their overlaps as small 
parameters; subsequently, it was shown by Read and 
Sachdev |l^ that the purely kinetic energy (T) piece de- 
scribed the 1/A^ dynamics of the nearest-neighbor SU{N) 
Heisenberg magnet in an extreme quantum limit. The 
former derivation is readily generalized to the triangular 
lattice and, crucially in this case, yields t > which we 
assume in the remainder; the latter is specific to bipartite 
lattices. 

T 3> V, t: At high temperatures, but less than the gap to 
non-valence bond states, static properties are obtained 
by a classical sum over all dimer configurations. Most 
crisply, consider T = oo where equal time correlators are 
given by unweighted averages. The square lattice prob- 
lem is critical in this limit, with algebraically decaying 
dimer-dimer correlations |T5[ , whence it is not surprising 
that at T = it orders everywhere except at a point. 

The important observation that motivated the present 
work, is that the triangular lattice problem is disordered 
in this limit with exponentially decaying dimer correla- 
tions. For example, the correlation function for two par- 
allel dimers separated by distance x along a column can 
be computed by standard Pfaffian/Grassman methods 
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FIG. 1. Phase diagram of the quantum dimer model on 
the triangular lattice. The nature of the ordered phases is 
indicated above the axis. 
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to see that the integrands are analytic in a finite interval 
about real values of the momenta, whence the two Green 
functions, G{A;B}i decay exponentially with x. 

This feature already yields a high temperature RVB 
phase, and makes a liquid phase far more likely at zero 
temperature, to which we now turn. 
Topological Sectors: An important ingredient of the 
analysis of the square lattice problem is the existence of 
two integer winding numbers, for periodic boundary con- 
ditions, that are conserved by any local dynamics in the 
dimer Hilbert space In that problem it is believed 
that the simplest dynamics (the analog of our T) is er- 
godic within each of the 0{L'^) sectors. On the triangular 
lattice, the situation appears to be quite different. Re- 
peating the analysis for the square case uncovers only 
four sectors corresponding to different combinations of 
even and odd windings which will be equivalent in the 
thermodynamic limit. With the minimal dynamics of T 
we have found that large classes of configurations can be 
connected to each other with no further constraint. The 
notable exceptions are the "staggered" configuration in 
Fig. Hand its symmetry related counterparts, which have 
no flippable plaquettes whatsoever. But even here a local 
four dimer rearrangement can be shown to allow a con- 
nection to other "generic" states. We conjecture there- 
fore that a general local dimer dynamics (which may re- 
quire as little as the inclusion of a four dimer move in 
addition to T) will be ergodic in each of the four topo- 
logical sectors and that the dynamics in our model is 
nearly so with the exception of the dynamically discon- 
nected non-flippable configurations. Note finally, that 
the Perron-Frobenius theorem implies that the ground 
state of H in each sector is nodeless. 
RVB Phase, Vc < v < t : As in the square lattice 
case, the triangular lattice dimer model has a "Rokhsar- 
Kivelson" (RK) point at t; = t at which the ground states 
are equal amplitude superpositions of all dimer coverings 
in a given sector. To see this, note that a lower bound for 
the ground state energy is obtained by considering each 



plaquette individually. A non-flippable plaquette is anni- 
hilated by H, whereas a flippable plaquette has a poten- 
tial energy of v and a kinetic energy of, at best, —t, which 
implies Eq > min{0, A'p(ti — t)}. The equal amplitude 



state \RK) = Nc J2c=i 1*^) sector) has energy 

{RK\ H \RK) — {v — t) (ufi) which vanishes and satu- 
rates the lower bound at v ^ t. Following our statements 
on the sectoral organization, we conclude that with the 
exception of the non-flippable configurations which triv- 
ially saturate this bound, there are four topologically de- 
generate ground states in the thermodynamic limit. The 
sum over all ground states is, for the purposes of comput- 
ing correlations diagonal in the dimer basis, equivalent to 
the classical dimer problem. As the staggered configura- 
tions are irrelevant to this sum, we conclude that the four 
generic sector states have exponentially decaying dimer 
correlations — i.e. they are RVB states. 

We next wish to argue that these states are repre- 
sentative of a phase at T = 0. We will argue shortly 
that their coexistence with the staggered states is due 
to a first order transition out of the RVB phase exactly 
at the RK point — much as in the square lattice prob- 
lem. Consequently, we wish to establish that there is 
a range Vc < v < t over which the RVB character of 
the ground state persists. Typically, one might antici- 
pate that a disordered ground state goes along with a 
gap to local excitations (as opposed to the degeneracy 
with globally distinct states). We have examined candi- 
dates for collective modes in the single mode approxima- 
tion and found that they are all gapped. We note that, 
in contrast, Rokhsar and Kivelson found that their criti- 
cal RVB state supported gapless excitations they dubbed 
resonons |17|. These together — the disordered character 
of the ground state and the gap in the local excitation 
spectrum — are strong evidence that the RK point is part 
of an RVB phase which is displaced to its right but will 
persist to its left (in Fig. |^). 

To test this argument, we have carried out quan- 
tum Monte Carlo simulations on systems upto 36 x 36 
sites at temperatures as low as t/10 using the method 
described in Ref. |^. As expected, we find that the 
dimer correlations are very short ranged and practically 
those of the classical dimer problem, conservatively, for 
2/3 < v/t < 1. We also find a very weak temperature 
dependence, suggestive of a gap. As v is reduced further, 
the correlations begin to exhibit crystalline structure (see 
below) . At V = t we find a hysteretic (first order) transi- 
tion to the staggered phase described next. 
Staggered Phase, v > t: For v > t the lower bound de- 
rived previously implies Eq > 0. As the staggered states 
are zero energy eigenstates of H they are the ground 
states in this range. This is similar to what happens in 
the square lattice problem, but the degeneracy is much 
lower in our problem (0(L°) vs 0{L)) and the exactness 
of the states is here a consequence of dynamics rather 
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than topology as mentioned earlier. 

An interesting consequence of this last observation is 
that the excitation gap is 0{L''^) in system size for the 
staggered phase At u/t = oo, the lowest energy ex- 
citations are the four dimer loop rearrangement of Fig. |^. 
At finite but large v/t, these get dressed by additional 
plaquette flips and acquire a dispersion. This branch 
of solitons is expected, by continuity, to be the relevant 
set of low energy excitations in the staggered phase even 
close to the RK point. 

Columnar Phase, — v t : To complete the phase di- 
agram, we move leftwards from the RVB phase in Fig. |^, 
turning first to the extreme case v/t — — oo. In this limit, 
the kinetic term T is disregarded, and the ground states 
are the maximally flippable states, i.e. those states |c) 
with maximal n/;(c). To identify the maximally flippable 
states, we note that each dimer can be part of at most 
two flippable pairs. Since the total number of dimers is 
fixed, the maximally flippable states are those in which 
each dimer belongs to two pairs. 

All maximally flippable states, whose number is expo- 
nential in L, can be obtained by carrying out the two 
operations A and B on a particular maximally flippable 
state as depicted in Fig. ||. Pairs of operations of either 
type as well as global translations and rotations can be 
shown to be generated by local plaquette flips. This illus- 
trates the point made earlier about winding number sec- 
tors. All maximally flippable states differing by an even 
number of operations are in the same sector, whereas a 
single A or B operation generates a state in another. 

Turning to the case of large but finite —v/t, we con- 
struct a perturbation theory using the small parameter 
t/v. Since any two maximally flippable states differ in at 
least 0{L) dimers, the degenerate perturbation theory is 
diagonal at any finite order. The energy shift of a state 
|c) at order 2n in perturbation theory depends on the 
number of flippable plaquettes, nfi{c'), of the states |c') 
which can be reached from |c) by at most n plaquette 
flips. 

The result of this perturbation theory is a strik- 
ing example of the phenomenon of quantum "order by 
disorder" — we find that it selects the columnar state de- 
picted in Fig. ||. States obtained by operations of type A 
are disfavored at 4"^ order in perturbation theory, those 
generated by type B operations at 6'^ order. 
Transverse Field Ising Point, v = : Together with 
P. Chandra we have recently shown that there ex- 
ists an exact correspondence between the quantum dimer 
model a.t V — 0, and the fully frustrated transverse field 
Ising model (FFTFIM) on the dual hexagonal lattice at 
fields r much smaller than the magnitude of the exchange 
J | [T9[ |. In this limit, the quantum ground state is con- 
structed entirely out of the ground states of the classi- 
cal frustrated model and the latter are (upto Ising de- 
generacy), in unique correspondence with dimer cover- 
ings of its dual, triangular, lattice. In our analysis we 




FIG. 2. Left: The columnar dimer state. The elementary 
dimer plaquette move generated by T is indicated in the bot- 
tom left plaquette. Such plaquette moves conserve the differ- 
ence between the number of dimers in rows marked by pluses 
and minuses. Dimer moves A and B, consist of shifting dimers 
onto the fat and dot-dashed bonds, respectively. Right: The 
staggered state with the four-dimer move connecting it to 
other states. 



found a low temperature crystalline "vl2 x vl2" phase 
which exhibits, in dimer language, a triangular super- 
lattice with a 12 site unit cell, that is consistent with 
a Landau-Ginzburg analysis of the Ising model. While 
the ordering observed at accessible system sizes is not 
conclusive regarding the fate of the model at T = 0, it 
appears that the columnar phase gives way to a different 
crystalline phase in the proximity of v ~ 0. 
Spinons: As noted before, the RVB phase has a gap to 
collective excitations. That is also true of the crystalline 
phases. Also of interest is the question of confinement 
for spinons — the gapped spin 1/2 excitations produced 
by breaking a valence bond. In the dimer model, these 
are represented by monomers or holons that carry a spin 
and so the questions of holon and spinon confinement are 
identical. This question is most easily addressed by con- 
sidering the free energies of states in which two monomers 
are held a fixed distance apart. At high temperatures, 
this is again a classical computation and it is clear that 
the spinons are deconfined on the triangular lattice . 
The contrast with the square lattice is again instructive, 
for there the spinons are confined at high temperatures. 
(However the confinement is very weak, only logarithmic 
|15| , which explains why the more disordered triangu- 
lar lattice does not confine.) At T = one can readily 
show that the state with an equal amplitude sum over 
dimer configurations with two spinons localized a fixed 
distance apart is an eigenstate at the RK point with an 
energy independent of their separation — i.e. the spinons 
do not interact beyond one lattice constant . By con- 
tinuity, we expect spinons to be deconfined in the entire 
RVB phase at T — 0. From these considerations it also 
follows that the charged excitation created by removing 
an electron from the system, will decay into a spinon and 
holon. Evidently, the holons/spinons will be confined in 
the crystalline phases. 

Spinon Confinement Transition: Several authors 
have suggested that a spinon confinement-deconfincmcnt 
transition will be governed by an Ising {Z2) gauge theory 
pO| , p^ 13 1 . In our own work we have found that frus- 
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trated transverse field Ising models (whose duals are pre- 
cisely the Ising gauge theories) can provide a description 
of valence bond phases of Heisenberg antiferromagnets on 
their dual lattices via their connection to quantum dimer 
models, such as the one considered in this paper. In the 
current context, the geometrical identification used by us 
leads to an intriguing observation. As already noted, the 
quantum dimer model is equivalent to the T <^ J limit of 
the FFTFIM at w = 0. Interestingly, the paramagnetic 
ground state of the FFTFIM at F ^ J, when projected 
onto the dimer manifold^ is the equal amplitude RVB 
sum that is the dimer model ground state at v = t. This 
suggests the conjecture that passing between these lim- 
its in the (projected) FFTFIM gives a description of the 
spinon confinement and translational symmetry breaking 
transition that takes place at the boundary of the RVB 
and \/T2 x \/T2 phases. We have argued in Ref. that 
the transition in the unprojected model is in the 0(4) 
universality class. It remains to be seen whether this 
identification survives projection. (The remaining tran- 
sition, between the columnar and \/T2 x \/T2 states must 
be first order on symmetry grounds.) 

In closing, we note that there are two previous "sight- 
ings" of a spin liquid phase on the triangular lattice in 
the literature. First, there is the large N, Sp{N) analy- 
sis of Sachdev |Q which found a disordered ground state 
with gapped unconfined spinous at small "spin" . Second, 
there is the exact diagonalization work of Misguich et al. 
p3[ on an 5 = 1 /2 system with a ferromagnetic two spin 
exchange frustrated by an antiferromagnetic four spin ex- 
change. They found a four-fold disordered ground state 
with a spin gap, but possibly confined spinous. We hope 
to clarify the connection between our results and these 
in the near future. Finally, we note that along the lines 
of the analysis in Q we expect doping to give rise, via 
holon condensation, to a superconducting phase on the 
triangular lattice. 

As we were finishing this paper, there appeared 
Rcf. lU], which reports neutron scattering evidence for 
deconfined spinous on an anisotropic triangular lattice. 
We note that the classical dimer problem in that case is 
also disordered and that one can construct an RK point 
with anisotropic terms that exhibits a disordered ground 
state. This suggests that the results of |24j could be un- 
derstood along the lines of our analysis in this paper. 
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